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SUMMARY 


A generalized linear surface noise interaction problem 
is formulated. Noise production by an oscillating surface, 
turbulent or vortical interaction with a surface and scattering 
of sound by a surface are Included in the generalized treat- 
ment. The direct effect of viscosity is included in the 
perturbation equations and the boundary conditions at the sur- 
face. An energy equation is derived to Illustrate how the 
work or virtual work done at the interface is partitioned 
into the acoustic and vortical modes. 

The problem is re-expressed in terms of a near field 
Incompressible viscous problem for the Bernoulli enthalpy and 
a problem for the compressible but Inviscid acoustic potential. 

A viscous Integral equation for the Bernoulli enthalpy, (or 
surface load) is derived. It is shown that viscosity will 
lead to a "unique” solution of the integral equation without 
the Kutta condition or other auxiliary singularity require- 
ment. The acoustic enthalpy is calculated by quadrature 
over the near field solution. 

The results of inviscid two-dimensional airfoil theory 
are used to discuss the interactive noise problem in the 
limit of high reduced frequency and small Helmholtz number. 

The acoustic spectrum is directly proportional to the surface 
load. The extreme limits of a full Kutta condition and no Kutta 
condition (actually no vorticity production) are considered. 

It is shown that in the case of vortex interaction with the sur- 
face the noise produced with the full Kutta condition is 3 
less than the no Kutta condition result. Also, the spectrum with 
Kutta condition decays monotonlcally with frequency while the 
corresponding result without the Kutta condition decays in an 
oscillatory manner. It is suggested that the difference in the 
high frequency spectra could be detected experimentally. 

The results of a supplementary study of an oscillatory 
airfoil in a medium at rest are discussed. It is concluded 
that viscosity can be a controlling factor in analyses and ex- 
periments of surface noise interaction phenomena. It is fur- 
ther concluded that the effect of edge bluntness and viscosity 
must be included in the problem formulation to correctly calcu- 
late the interactive noise. 



I. INTRODUCTION 


In two previous reports (Refs. 1 and 2) the homen- 
troplc theory of aeroacoustlcs , based on a kinematic 
definition of sound and the concept of Bernoulli enthalpy, 
was formulated in detail and applied to a variety of pro- 
blems that illustrate the three basic questions of ’’free 
flow" aeroacoustlcs (i.e., in the absence of solid boun- 
daries). In the spirit of Chu and Kovasznay (Ref. 3) the 
modal concept has been stressed in all of the previous 
work. The problem of how to decouple the acoustic and 
vortical modes (to the extent that it is possible) of 
energy transport is one of the most interesting questions in 
the field of aeroacoustlcs. 

For example, when the problem of noise production by a 
turbulent flow (e.g., a jet) is considered, it is customary 
to adopt the "Lighthlll hypothesis". Because of differences 
in the basic definition of sound in the various aeroacoustic 
formulations and the level at which the Lighthill hypothesis 
is invoked, there have been serious arguments over what con- 
stitutes the "source" of sound. With the modal approach 
there is no "source" per se, and the application of the Light- 
hill hypothesis is tantamount to assuming a unilateral trans- 
fer of energy from the vortical mode to the acoustic mode. 

The only internal mechanism for acoustic feedback to the vor- 
tical mode is the corlolis coupling acceleration between the 
acoustic particle velocity and the vorticity of the primary 
flow (see Eq. (2.28)of Ref. 1). The corlolis mechanism is 
Important in the problem of sound scattering by a steady vor- 
tical flow and the stimulated emission of sound from an un- 
steady vortical flow. Both of these problems were inves- 
tigated in detail in Ref. 2. It was shown by explicit calcu- 
lation that the broadband noise radiated from a discrete 
vortex array can be enhanced by an incident sound field. 

(For a concise summary of this work, see Ref. ^). 

The distinction between acoustic and vortical energy 
transport and the corlolis coupling mechanism between the two 
modes is reasonably well understood in the description of the 
internal dynamics of a fluid medium. However, there remains 
the important problem of how these modes Interact in the 
presence of solid boundaries, and in particular when the 
solid boundary has "sharp" edges. The Important feature of 
the surface interaction problem is that viscosity becomes a 
controlling factor in the mode coupling at the surface of 
the body. It is the objective of this report to shed some 
light on the importance of viscosity in a class of aeroa- 
coustic problems that involves thin airfoil like surfaces, 
with edges. Two experiments due to Brooks (Ref. 5) and Brooks 
and Hodgson (Ref. 6) will be used to guide and test the 
development of the theory. 
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NOMENCLATURE 


a 


A 


n 


c 

c 


« 


C 

mn 

E 

f 


F 

G 

h 


h ' 




Jn(z) 

k 

k 


Ko(z) 

1 


M 


free stream speed of sound 

spectral coefficients of the surface load, 
see Eq. (2.8l) 

airfoil chord 

compressibility, see Eqs. (2.72) and (2.73) 

denotes chord of the surface S in the two- 
dimensional problem 

see Eq. (2.83) 

energy density, see Eq. (2.9) 
surface deflection along z-axis 
surface force distribution, see Eq. (2.47) 
magnitude of the surface force, see Eq. (2.47) 
see Eq. (2.84) 

Green’s function, see Eqs. (2.75) and (2.7^) 

vertical distance of vortex from interactive 
surface 

perturbation enthalpy 
Bernoulli enthalpy 

unit vectors along x,y,z axes, see Fig. 1 

Bessel function 

u)c/2Vj^, reduced frequency 

oj/ a 

modified Bessel function 

spanwlse correlation length of the surface load 
Vop/a, free stream Mach number 


-*■ 

n 

q 

Aq 


Q 

R 

S 

S(a) 

t 

Tn(x) 

Un(x) 

V 

oo 

V 
W 

x,y ,z 


raomentura flux, see Eq, C2.10) 
unit normal to a surface 

denotes Fourier transform of variable q with 
respect to t5.me 

q(z = 0") - q(z = O'*'), Jump in the dependent 
variable q across the plane z = 0, e.g., see 
Eq. (2.53) 

dissipation function, see Eq. (2.11) 
see Eq . (2.77) 

denotes a fixed or oscillating surface 
outer boundary of the volume V 

acoustic spectral function, see Eqs. (3-1) and 
(3.2) 

time 

Chebyshev polynomial of the first kind 
incompressible perturbation velocity 
Chebyshev polynomial of the second kind 
perturbation velocity field 
disturbance velocity field 
free stream velocity 

integration volume in energy integral 

upwash function see Eq. (2.^) 

Cartesian coordinates, see Fig. 1 

acoustic intensity, see Eqs. (3-36) and (3*37) 

kernel function, see Eqs. (2.65) and (2.66) 

integrated kernel, see Eq. (.2.85) 

total lift on two-dimensional surface, see Eqs. 
( 3 . 16 ) and (3.10) 

normalized surface load distribution 
surface work or virtural work, see Eq. (2.13) 
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IT 

o 

a 

y(x) 

r 

6(z) 

e 

K 

x,x* 

A 


V 


T 

<P 

X 

u 


curl 

div 

grad 

Vi 


see Eqs . (3*28) and (3-3^) 
see Eq. (3-3) 

1 - 

vortex strength distribution, see Eqs. (3.9) and ( 
circulation, see Eq. (3.16) 
vortex strength, see Eq. (3-33) 

Dirac delta function 

far field directivity angle 

o)c/2a, Helmholtz number 

see Eq. (2.57) 

see Eq. (2.79) 

kinematic viscosity 

see Eq. (3.35) 

acoustic potential 

hydrodynamic potential, see Eq. (3-7) 
frequency in Fourier transform, see Eq . (2.^2) 
perturbation vorticity, see Eq . (2.21) 

vector curl operator 
divergence 

vector gradient operator 
Laplace operator 

surface Laplace operator, see Eq. (2.60) 
time average 
absolute value 


3.15) 
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II. GENERAL THEORY OF A CLASS 
OF LINEAR AEROACOUSTIC PROBLEMS 

A. The Surface Interaction Problem 

Below, a linear problem is formulated that will permit 
the investigation of three basic noise-surface interactions: 

1) An oscillating surface; 

2) An assumed vortical disturbance (e.g., turbulence) 
interacts with the surface (The Edge Noise Problem); 
and 

3) An assumed acoustic disturbance interacts with the 
surface (The Diffraction Problem). 

These problems are depicted schematically in Figure 1. The 
surface S may be finite or infinite and initially the sur- 
face is supposed to be of zero thickness. Later, in Section 
III, the effect of finite edge geometry, is discussed. 

It is assumed that each of the three problems may be 
described with the theory of homentroplc small disturbances 
from a uniform main stream. The appropriate perturbation 
equations are: 

^ Dhl ^ dlv7- = 0 (2.1) 


+ gradh' = vV^v’ + ^ grad (divv') 

where 

A = _1_ + V -L 
Dt 3t “ 3x 


( 2 . 2 ) 

(2.3) 


For the problem of an oscillating surface the boundary condi- 
tions on S are as follows: 


V' = W = ^ 


Df 

Dt 


on z = 0- in S 


k X v’ - 0 


( 2 .^) 


where f is the transverse (along the z-axis) deflection of 
the surface. The boundary conditions for the other two problems 
can also be expressed in the form of Eq. (2.4). Denote the 
incident velocity field due to ^ vortical or acoustic distur- 
bance in the absence of S by Vj. Then the interactive 
disturbance field must be such that it cancels the velocity v. 
on S; i . e . , ^ 
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Noise Produced 



Figure 1 -> Three Types of Interactive Noise 
Production Mechanisms 
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-> -4- -V 

K»v’=W--k-v^ 
k X V* = “k X 


on z 


0± In S 

(2.5) 


The boundary condition on the vertical, or z-component 
of velocity, is the usual boundary condition that is familiar 
in inviscid aerodynamic problems. The second boundary condi- 
tion on the tangential surface perturbation velocity is necessary 
when viscous terms are retained in the momentum equation. The 
importance of the viscous terms and boundary conditions in 
establishing uniqueness of solution (without a Kutta condition) 
and in the partitioning of disturbance energy into the acous- 
tic and vortical modes is the essential content of this report. 
Because of the linearity of the boundary value problem the 
solution can be obtained by superposition of two parts with 
the respective boundary conditions: 

k • v’ = W = -k • v^ 

ic X v’ = 0 

and 

"I 

k • V ' = 0 

on z = 0- in S 

5xv’=-icxv^ (2.7) 


on z = 0- in S 

( 2 . 6 ) 


For all finite disturbance velocity fields v^ the solution 
of Eqs. (2.1) and (2.2) with the boundary conditions Eq. (2.7) 
is small of the order of the viscosity and will not be con- 
sidered further in this report. The homogeneous viscous 
boundary condition in Eq. (2.^) or (2.6) is the main result 
that is needed to formulate The generalized linear aero- 
acoustic problem. Solve the perturbation equations (2.1) 
and (2.2) subject to the boundary conditions (2.4) or (2.6) 
and the condition of outgoing decaying disturbances at Infinity. 


B. A Self Consistent Energy Equation 

It is instructive to derive an energy equation for second 
order perturbation quantities. In so doing the global par- 
titioning of energy into acoustic and vortical parts will 
become evident. Multiply gq. (2.1) by h’, form the vector dot 
product of Eq. (2.2) with v' and add the results to obtain the 
local energy balance. 


DE 

Dt 


+ dlv M 


- -Q 


( 2 . 8 ) 
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where 


E = + 

v’ ^ 

(2.9) 

2a^ 

2 

1*1 = h’v’ -vv* X co’ 

11 ^ 

- -jvv ’ div V ’ 

(2.10) 

Q = v|w’|^ + 

(div V ’ ) ^ 

(2.11) 


For a statistically stationary disturbance velocity field, 
the time average (denoted by an overbar) of Eq . (2.8) is 

Voo|| + div5^ = -Q (2.12) 

Now integrate the last result over a large spherical volume, 

V, that encloses the surface S. On the outer boundary, S^, 
the effect of viscous transport is assumed to be negligible. 

The resultant global energy balance becomes 




+ dxdy 




h ' dA + 


Voo? * n E dA 
QdV (2.13) 


For the case of an oscillating airfoil, the first example of 
the generalized problem, the quantity ^^is the rate at which 
work is done on the fluid by the oscillating surface. For 
the second and third examples the surface is rigid and does 
no work on the fluid. However, it is convenient to interpret 
^^as the virtual work due to the surface motion (upwash) W 
that must be imposed to counter the normal component of the 
disturbance velocity field. The global energy equation (2.13) 
shows that the work done by the surface S is partitioned 
into sound radiation through the outer surface S^, convective 
energy transport through Sf, and viscous dissipation through- 
out the volume V. The last two parts are actually different 
manifestations of the same physical process; i.e., the forma- 
tion of vorticlty at the surface S. If the surface Sf is 
sufficiently far from S, no convective transport will be de- 
tected and all of the vortical energy will be counted in the 
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viscous dissipation. If is closer to S then convective 
transport of vortical energy will be a significant part of 
the energy balance. It is noted that similar statements can 
be made about the acoustic term in Eq. (2.13). Certainly 
if Sf is removed to infinity then all of the surface work is 
dissipated by viscosity. 

Finally it is pointed out that viscous dissipation cannot 
be discarded as an energy sink by the simple argument that 
the viscosity is small. It will be shown by explicit cal- 
culation in Section III and in the supplementary report that 
the vorticity is sufficiently singular to yield a finite in- 
tegrated value for the dissipation even for vanishing small 
viscosity . 

C. Mode Splitting and Boundary Conditions 

The boundary value problem composed of equations (2.1), 
(2.2) and the generalized boundary conditions (2.^) or (2.6) 
is well posed. It is of interest, however, to recast the 
problem in terms of acoustic and vortical modes as follows: 

Let 


v' = grad (J) + u’ (2.1^) 

where (p is the acoustic potential and u’ is an incompressible 
velocity field. Substitute Eq. (2.l4) into (2.2) to obtain 

^ + grad'H = (2.15) 

where 


is the Bernoulli enthalpy (or Pseudo Sound) associated with 
the incompressible problem for 3'. Now use Eq. (2.16) to 
eliminate h^ in Eq. (2.1). The result is a wave equation for 
<p f i . e . , 


a2 Dt2 


4 

3 a2 


7254 - 

Dt Dt 


(2.17) 


Note that, equations (2.15) and (2.17) could also be obtained 
by direct linearization of Eqs. (2.27) and (2.28) in Ref. 1. 
The acoustic mode is "driven" in the near field by the 
substantive rate of change of H . In the far field the 
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acoustic mode decays due to viscous dilltational damping. 

The near field hydrodynamic or vortical mode is not coupled 
linearly to the acoustic mode In the interior of the gas. 
However, mode coupling does occur due to the boundary condi- 
tions at the surface. (See Ref. 7 for further discussion of 
this point). 

The boundary conditions on u'jiy are chosen to be the 
same as those for the complete perturbation velocity field; 

1 . e . , 


k • u ’ = W 


on z = 0- in S (2.l8) 
Then the boundary conditions for <p are homogeneous; i.e.. 


Ic X u’ = 0 


k • grad(f) = 0 
k X grad(J) = 0 


on z = 0- in S (2.19) 


It is possible to apply the no-slip boundary conditions on 
3’ because of the viscous term in Eq. (2.15). Similarly it 
is possible to apply no-slip conditions on the acoustic par- 
ticle velocity as long as the viscous damping term is retained 
in Eq. (2.17). It is the presence of these terms that pre- 
cludes the necessity for a Kutta condition or other singularity 
criteria to solve the relevant boundary value problem. 


A somewhat different insight into the role of viscosity 
and the coupling of the pressure and vorticity modes may be 
obtained as follows. Take the curl of Eq. (2.2) or (2.15) to 
obtain the vorticity diffusion equation 


Dw’ 

Dt 




( 2 . 20 ) 


with 


= curl v * = curl u 


( 2 . 21 ) 


Also take the substantive derivative of Eq. (2.1) and the di- 
vergence of (2.2) to obtain the damped wave equation for h’; 

i.e. , 


11 


3 a2 Dt 


( 2 . 22 ) 


1 

a2 Dt2 


- V^h» *= 0 


As noted by Chu and Kbvasznay (Ref. 3) the vorticity and pres- 
sure modes decouple completely in the interior of the gas. 

At the surface, the boundary conditions on can be used to 
derive the relations between h’ and 3'. First write Eq. (2.2) 
in the form, 

W + srad(h' + i ^ ^ - vcurlS- (2.23) 

Because of the planar geometry of S and the no-slip boun- 
dary conditions, it follows that 

032 = ic • = 0 (2.24) 

Now project Eq. (2.23) onto the surface S and use the 
boundary conditions on to obtain the following boundary 
conditions . 


4 

V 

Dh' 

\ - DW 

3 

a^ 

Dt , 

1 " Dt 

4 

V 

Dh’ 


3 

a2 

Dt , 

) Bz 


9o)^ 9g)y 
By ” Bx 


on z = 0“ in S 


(2.25) 


The intrinsic coupling of the pressure and vorticity modes 
is evident from these relations. The normal pressure gra- 
dient at the surface is balanced by the surface acceleration 
(real or virtual) and gradients of the surface vorticity. 

The second expression shows how vorticity is produced at a 
surface by pressure gradients along the surface. Near sharp 
edges that have large pressure gradients, the production of 
vorticity can be Intense. It is also clear from this rela- 
tion that to obtain a meaningful solution (non-constant 
surface load) of the viscous boundary value problem, the pro- 
duct v( Sw’/Bz) must be finite even for a vanishingly small 

viscosity coefficient. The coupling between modes is thus 
very singular. 
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D. The Associated Inviscld Boundary Value Problem 

For each of the boundary value problems and energy rela- 
tions derived thus far there is a corresponding "inviscld" 
relation. First of all, the perturbation enthalpy and 
velocity satisfy the relations 


a 2 Dt 


+ dlvv’ 


0 


( 2 . 26 ) 


™ + gradh' = 0 (2.27) 

with the boundary condition 

ic • v' = on z = 0^ in S (2.28) 

and the radiation condition. The energy relation (2.13) 
becomes 





n "E dA 
(2.29) 


where the outer surface can be taken to infinity since 
acoustic and vortical waves do not decay. The convective 
transport of energy in Eq. (2.29) is due to the singular vor- 
tex wake. It is the singular remnant of the effect of vis- 
cosity, and in general it is a significant part of the work 
done by the surface. Without the explicit effect of viscosity 
at the surface there is no way to determine the energy parti- 
tioning without recourse to an artificial uniqueness criteria 
like the Kutta condition. This will be demonstrated explicitly 
in the next section where the boundary value problem is cast 
in the form of a singular integral equation. 

The inviscld boundary value problems for u’ and (|) are 
respectively : 


div u' = 0 


(2.30) 


Du’ 

Dt 


+ grad'H = 0 


(2.31) 
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( 2 . 32 ) 


and 


ic • u’ « W on z = oi in S 


with 


Rli 

a2 Dt2 


- V^<{> = 


a 2 Dt 


(2.33) 


li 

8z 


0 on z = oi in S 


(2.3^) 


The perturbation enthalpy 


h' 


'H - 


D4> 

Dt 


also satisfies the homogeneous wave equation 


(2.35) 


1 D^h’ 
a2 ^t^ 


V^h' = 0 


(2.36) 


and 


- 4 - 

Deo' 

Dt 


= 0 


(2.37) 


It follows from (2.25) that 


™ on z = 0± in S (2.38) 


In each boundary value problem the boundary conditions on 
the tangential surface velocity components must be given up 
because the inviscid equations are of lower order-. Also, the 
inviscid problems do not have unique solutions and some 
auxiliary criterion (usually the Kutta condition) must be 
used to establish uniqueness. 

E. An Integral Equation for Bernoulli Enthalpy (Uniqueness) 

Consider the viscous boundary value problem for u’ and 
Ti ; l.e. , 
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dlv u’ 


0 


(2.39) 


Du' 

Dt 


+ grad"H 


vV^u' 


(2.40) 


ic • u' = W 
X u' = 0 


on z = oi in S 


(2.41) 


It Is convenient to work in the frequency domain. Take 
Fourier transforms of the above equations with respect to 
time and note for any dependent variable, q, that 


q(o)) = q(t)dt 

q(t) = i qCo) dii (2.1|2) 

Then 


div u* = 0 


(2.43) 


with 


and 


+ grad;K = vV^u' 


(2.44) 


S . u' = W 

ic X u' = 0 


+ 

on z = 0- in S 


(2.45) 


D 

Dx 


3 

3x 


+ io) 


(2.46) 


15 


The solution can be represented in terms of an unknown nor- 
mal force distribution on the surface S; l.e., 

I =• -S FCx,y) 6 (z) ( 2 .'I 7 ) 

Add ^ to the right hand side of (2.AA) and then take the 
divergence and curl to obtain 

= -F -P. (2.48) 

~ ^ dz 

and for the vorticity 

Dw’ ^ 

w’ = k X gradF6(z) (2.49) 

The solution of (2.48) is 


■V = 


1 9 

¥tt 3z 



Z(?) d? 
lx - ll 


where the integration is over the surface S with 


( 2 . 50 ) 


f = Jc + Jn , d| = dc dn (2.51) 

Take the limit of (2.50) as z 0 to obtain 


+ F 

2" = Limit =+^ (2.52) 

z -► 0- 


or 


F = Alt = -lj(x,y,0+) (2.53) 

is the local jump in the Bernoulli enthalpy on the surface. 
For points off the surface must vanish. 

The solution of (2.49) can be expressed in the form 

3’ = £ X grad Q (2.54) 
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where 


(2.55) 


W -V V^Q = F 6(z) 


and It Is readily shown that 


Q = 


_ 1 


TtFv 


f p(f) exp[xt • (^ - I) - x»|; - g|] 


5 - 5l 


where 


and 


= (X^ + , X = 


Voo 

2v 


curlw’ = itV^Q - grad 

o Z 


From (2.^4), with z ^ 0 


dC 

(2.56) 

(2.57) 

( 2 . 58 ) 


where 


Also 


Thus 


DW , 9:H ^ 

= -vk • curlu) ’ 

Dx 3z 

= -vVjQ 


2 _ 


9x2 9y2 


12 . 

9z 


= ’1 /s 


F(C) 


X --Cl 


dC 


(2.59) 


( 2 . 60 ) 


( 2 . 61 ) 


DW _ 1 

Dx TTif 


/s 


1 - exp[x? • (x - C) - X*|x - fl] 


dt 


X - C 


for z 0 (2.62) 

Now solve (2.62) for W and take the limit as z 0 to obtain 
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(2.63) 


the following integral equation for f(C)j i.e., 

£(t)JT(x - f)at = W(|) , X in S 

Where the kernel is defined by the relation 


^(X) 



V 


2 

1 



exp(-lo)s/v^) 


1 - exp{X(x - s) - X*| 


X - 


Ts 



ds 


C2.64) 


The integral equation for the corresponding two-dimensional 
case can be obtained by integrating (2.64) over the spanwise 
coordinate. From a previous report (Ref. 8) the final ex- 
pression is 

y £(5)^(X - 5) dc = - w(£) . X in C (2.65) 

Where 

Tlnlx - £| + exp{X(x - £)} K^^(X*|x - £|)jd£ 

( 2 . 66 ) 

and the integration is only over the chord C of the inter- 
active surface S. 

The solution of (2.65) and a more general equation for 
airfoil type surfaces with finite thickness has been the sub- 
ject of a recent development of "viscous thin airfoil theory". 
The results are summarized in a final report (Ref. 8) and an 
AGARD report (Ref. 9)* For the moment the most important 
point to note is that the kernel function ,_^^(x) (Eq. (2.66)) 
has a logarithmic singularity for small argument no matter how 
small the coefficient of viscosity. The corresponding Invls- 
cid kernel of the incompressible Posslo equation (obtained 
by dropping the second term in (2.66) has a Cauchy singularity. 
The viscous integral equation has a unique solution (See 
Ref. 8) while the inviscid equation has an eigensolutlon that 
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Is associated with the circulation of the perturbation flow 
field (See Section IIIB) . The strength of the eigensolution 
determines uniquely the amount of vortical energy that is 
shed into the wake. The eigensolution also contributes di- 
rectly to the "source" term Dt/ /Dt in the acoustic problem 
( 2 . 33 ) and ( 2 . 3 ^) and so affects the radiated sound. With- 
out viscosity there is no physical criterion to determine 
the eigensolution and it is customary to resort to the Kutta 
condition or principle of minimum singularity at the trail- 
ing edge to establish uniqueness. The solution of the 
viscous Integral equation (2.66) has a square root singular- 
ity at a mathematically sharp leading or trailing edge. The 
singularities can be removed by slightly blunting the edges 
as shown in Section III and in that case there is no useful 
edge criteria whatsoever to establish uniqueness. The sur- 
face load, shed vorticity and noise must be obtained by 
solving the viscous boundary value problem. 


F. The Acoustic Potential 


Once is known from the solution of the appropriate 
Integral equation, there remains the problem of solving 
for the acoustic potential l.e.. 


with 


-L 

a2 Dt2 


ii 

3 a2 Dt 


V^4) 


1 

a 2 Dt 


(2.67) 


3 (j) 

Jz 


= 0 


k X gradcf) = 0 


on z 


0- in S 


( 2 . 68 ) 


and the radiation condition. Actually, because of the dila- 
tational damping term, it is sufficient to require decaying 
solutions of ( 2 . 67 ) at infinity. 


The Interpretation of (2.67) and (2.68) is conceptually 
appealing. The origin of the interactive noise is the local 
substantive rate of change of V* The surface itself is a 
passive element in the noise production process but it does 
cause compressible diffraction primarily through the boundary 
condition on the Induced vertical component of velocity. Also, 
it is emphasized that by construction no vorticity is pro- 
duced in the compressible part of the Interaction process. 

All vorticity is associated with the surface loading due to 
'h. It is this observation that permits an important simpli- 
fication of the acoustic problem; l.e., the direct effect of 
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viscosity can be omitted in the acoustic problem. Since 
there must be no compressible production of vorticlty, the 
compressible viscous term is omitted in (2.67) and the viscous 
boundary condition in (2.68) is replaced by the condition that 
the acoustic potential be continuous across the plane z = 0 
for all points not in S; l.e.. 


a2 Dt2 


- = 


_L ^ 

a2 Dt 


(2.69) 


~ =0 on z = 0- for (x,y) InS (2.70) 

4>( z = O'*') = Hz = 0") for (x,y) not in S (2.71) 


plus a radiation condition on the far field behavior of (J) . 

Now take the Laplaclan of (2.69) and introduce the com- 
pressibility variable c* such that 

V2<j) = c* (2.72) 

Then, in the frequency domain 


A. dV 

a2 Dx2 


V^c* 


a 2 Dx ~ 



a 2 Dx 3z 


(2.73) 


where (2.^8) has been used to replace . The solution for 


where 


c* ^ 1. A. 1. 
~ a2 Dx 9z 


1 D^G 
a^ Dx2 


Js GCx 

- V^G - «(x) 


C) dc 


(2. 7^0 

(2.75) 


or 


g _ exp(-lk/6^ ) (R - Mx) 
iittR 


( 2 . 76 ) 
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with 


R = + 3^ (y^ + )j ^ 

3 ^ = 1 - 
M = v„/a 
k = co/a 


(2.77) 


In the far field 


c* = J:- A G(x) f exp(ikA/3^) F(|) df (2.78) 
a ^ Dx 9 z J S ~ 

with 

A - (I - n (2.79) 

From the continuity equation (2.1) it follows that the far 
field acoustic enthalpy is- the compressible continuation of 
V via the relation 


/ exp(ikA/gM F(f) d| (2.80) 


Thus it appears that the far field pressure and the OASPL 
for example can be calculated from the incompressible 
problem without specifically solving for the acoustic poten- 
tial. The actual acoustic potential and velocity and there- 
fore the far field Intensity must Include the solution of 
( 2 . 72 ) subject to the diffraction boundary conditions (2.70) 
and (2.71). 

The essential effect of viscosity, then, is to establish 
"uniqueness" of solution of the surface load distribution 
in the integral equation (2.65) with the kernel function 
(2.66). Once the surface loading is determined, the shed 
vorticity and the acoustic far field are uniquely determined. 
The same statement can be made about the corresponding In- 
vlscid problem (See Section IID). However, the corresponding 
compressible or incompressible integral equation for the sur- 
face load has an eigensolution . The strength of the acoustic 
field and the vortex wake depend on the strength of the 
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eigensolution. Without a viscous integral equation to de- 
termine the surface reactive load there is no apriori phy- 
sical criterion to determine the eigensolution. It is then 
customary to invoke the Kutta condition, principle of mini- 
mum singularity or modifications thereof to establish 
uniqueness. It is important to note that edge singularities 
are completely admissable in the linearized viscous or In- 
viscld problem. Only by accounting for finite edge curvature 
can these singularities be eliminated. It is also true that 
the edge singularities pose no real problem in the calculation 
of the noise or vortex wake. It is much more important to 
calculate the surface loading correctly, both in magnitude 
and phase, than to worry about the appearance or absence of 
an edge singularity. These remarks, by the way, are indepen- 
dent of frequency. Even for very low frequency interactions 
where the noise and wake vortlcity scale with the total lift, 
it is known (Ref. 8) that the Kutta condition will usually 
overestimate the lift. 

G. Chebyshev Representation of the Surface Load 
and Acoustic Far Field 

Consider the class of quasi two-dimensional problems 
where the surface load is calculated from (2.70) and the far 
field from (2.79) with the effective spanwise length scale 
/ over which the load is correlated or non-zero. The solu- 
tion of the integral equation (2.65) can be expressed in a 
series of Chebyshev polynomials; i.e.. 


F(x) 


00 


E 

n=0 


Tj^(2x/o) 

[l - (2x/c)^]^ 


( 2 . 81 ) 


where c is the chord of the two-dimensional surface. The 
spectral coefficients A^ may be calculated from the finite 

matrix representation of the integral equation; i.e., 

N 

V An = 6m m=0,l,..,N (2.82) 

with 

p _ (m + 1) A 

Cmn ^ J ^ 

(2.83) 


f ^ 

cosnc|)dtt) / cos{ (m+1) e } de,_^^(cos 0 - cos (J)) 

Jo ^ 
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and 


Ujn(x) W(x) dx (2.84) 






^U) <iK 


(2.85) 


The numerical solution of (2.82) for the surface load distri- 
bution is the subject of Ref. 8. 

Now substitute the Chebyshev series (2.81) into (2.80) 
to obtain the acoustic far field in terms of the spectral 
coefficients of the surface load. The result is 




— . f 
dz ^ 



exp [(iK/B^)(x/R - M )]? 


00 


Z 

n=0 


Tn(2C/c) 

" (l - (2C/o)^) 


d? 


3 G 

3z 


£c 

2 


°0 - 

/ exp{(i</eM 

n=0 Jo 


(x/R - M)cos9} cos n9d0 


_ i K * i 

8 


(z/R) exp{ (-1 K ) (R - Mx)} 


An 1^ Jn^( </B^)(x/R - M)} (2.86) 

n=0 

with the Helmholtz number 


< = o)c/2a = kc/2 (2.87) 

For the range of compact frequencies ( < << 1) only the 
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first term In the series of (2.86) contributes to the far 
field and the magnitude of h^ is proportional to or the 

total interactive lift that develops on the surface. The 
value of Aq may be calculated uniquely with the viscous 
theory as shown in Refs. 8 and 9. 

H. Effect of Edge Bluntness 

From extensive calculations of the solution of (2.82) 
in Ref. 9 the following important result was found. For a 
surface with zero thickness the edges are mathematically 
sharp. As a result, the solution of (2.82) converges to the 
Inviscid solution with Kutta condition (See Section III) for 
any Reynolds number greater than about 500. The departure 
from inviscid behavior is of order l/(Re)^ . The reason 
is that by introducing a mathematically sharp edge or zero 
length scale the linear viscous effect is grossly overex- 
aggerated at the edges. By introducing a finite edge radius 
the viscous edge effect is diminished and viscous action over 
the entire surface becomes the controlling factor in establish- 
ing the surface loading. In all cases calculated (See Ref. 8) 
the surface loading is diminished by Introducing finite edge 
geometry. An example of the reduction in the steady state 
lift curve slope on a surface at angle of attack is shown in 
Figure 2. The lift on an elliptic section is about 15^ less 
than the lift on a Joukowskl airfoil with cusped trailing 
edge. !More extensive comparisons of unsteady results for 
surfaces with and without edge curvature may be found in Ref. 

8. For moderately low values of reduced frequency (cjc/2voo = 
order of a few tenths) the real and imaginary parts of Aq, A]_, 

the two leading spectral coefficients in the surface load 
series (2.80) for an elliptic section are substantially 
different from those of flat plate theory. The changes are 
in quantitative agreement with measurements of Davis and 
Malcolm (Ref. 10). 

To perform meaningful acoustic calculations it is neces- 
sary to solve the surface load problem for relatively large 
values of reduced frequency { (cuc/2v») = 0(1/M) and larger }. 
The numerical problems in evaluating the kernel functions 
with thickness and viscosity for high frequency have not been 
completely resolved and detailed results on this aspect of 
the problem will be reserved for a future study. 

III. APPLICATIONS OF "INVISCID” THEORY 

A. Zero Thickness Two-Dimensional Surface 

In the following development the family of solutions of the 
well known two-dimensional incompressible oscillating airfoil 
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Figure 2 - Degradation of the Steady State Lift Curve Slope (C, ) 
due to Trailing Edge Bluntness (e*); note the a 

airfoils sketched above the graph. 



problem (Ref. 11) is used to draw some general Inferences 
about the generation of sound for the generalized surface 
interaction problems. Following the work of Yates and 
Houbolt (Ref. 12) and recent work of Howe (Ref. 13) the ex~ 
treme limits of a full Kutta condition and no Kutta condition 
are used to illustrate why it is important to complete the 
development of the viscous theory as formulated in the pre- 
vious Sections. 

For the two-dimensional problem, the spectrum of the 
acoustic enthalpy (See Eq. (2.79)) can be expressed in the 
form 

= - H “t- S(a) (3.1) 

where G is the asymptotic form of the three-dimensional Green’s 
function ( 2 . 76 ) and i is some measure of the spanwise corre- 
lation of the surface disturbances. The airfoil semi-chord 
(c/2) is chosen as the unit of length and the spectral func- 
tion S(a) is expressed in the form 


S(a) = J exp(iax) .2^(x) dx (3-2) 

where 

a = (iK/6^)(cos0 - M) (3-3) 


and _^(x) is the normalized chordwise load distribution 
that must be calculated from the incompressible boundary 
value problem: 


div u = 0 


with 


Du 

Dx 


+ grad P 


0 


j • u = ony = 0, -l<x<l 


(3.4) 


D 

Dx 


_ 3 _ 

3x 


+ ik 


(3.5) 
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and 


^(a) = P(x,0-) - P(x,0+) 


(3.6) 


The solution of (3.^) ca-n be e:xpressed in terms of a hydro- 
dynamic potential \ and a vortex sheet distribution y(x) 
that are related to as follows : 


u = grad X for y 0 


U) 


3v 

3x 


3u 


3y 


= y (x) 6(y ) 


C3.7) 

(3.8) 


y(x) 


/"* , 3 Ax 

= c, dy = -^ 


(3.9) 


- 


DAx 

Dx 


( 3 . 10 ) 


with 


AX = X(x,0“) - X(x,0''') 


(3.11) 


The potential is a harmonic function with the representation 


X = 


L f 

!7T J - 


1 Y(5)tan-1 (x ' - c ' ) 


dC (3.12) 


For X > 1 (in the wake) the following relations hold; 



^(x) = 0 

(3. 

13) 

AX 

= - rexp{-lk(x - 1)} 

(3. 

IM) 

y(x) 

= - lkrexp{-ik(x - 1)} 

(3. 

.15) 


where T is the negative of the potential Jump across the 
trailing edge. Also, r is the circulation around the 
surface ; i . e . , 
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y*o- 

r= JP U*ds = -Ax 


/ 

y«o+ 

X=1 


( 3 . 16 ) 


x=l 


The value of f is unknown and in fact is completely 
indeterminate with Invlscid theory. 

An integral equation for y(x) Is obtained by applying 
the boundary condition (3*^) to the integral representa- 
tion (3.12). The result is 




= S(X) 


(3.17) 


with 


g(x) = -2 2^(x) 




exp { -ik(g - 1)) 


(C - X) 




where the slash through the integral denotes 
the Cauchy principle value. 


(3.18) 


The solution of (3.17) is 


Y(x) = 


/: 


7T ( L - X^ ) 

exp{-lk(C - 1)} 




2 


(1 - s^ ) 
s - x 


ds 


U - x) 


+2 




(1 - S^) 




dC 


^^(s) ds 


(3.19) 


-1 (s - x) 


a well known result (e.g.. See Ref. 11). 

The spectral function S(a) from (3.2) can be expressed 
in terms of y(x) with the formula 
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S(a) 


r e^“ + (I 



tax 

e 


Y (x ) dx 


( 3 . 20 ) 

that is easily derived from (3.9) and (3.10) by integration 
by parts. Now substitute (3.19) into (3.20) to obtain the 
following representation of the acoustic spectrum: 


S(a) = r 


1 _ !S) j (a) + k glc 
a / o ^ ' a 


+ 21k 


- i" J„(a) f" exp{-lkt£ - 1)1 dg' 

“/ " y 1 [£+(£'- 1 


+ H 


CO p 

-t) § i" j\ a - 5^)’^ V,(C)^C?)d£ 


(3.21) 


where Jj^(a) is the ordinary Bessel function and Uj^(C) is the 
Chebyshev polynomial of the second kind. 

On the basis of the viscous analysis of a zero thickness 
surface, (Ref. 9) the value of the circulation T must be 
such that the Inviscid flow at the trailing edge is smooth 
(the Kutta condition). In that case the expression in square 
brackets in G.19) must vanish as x 1”. The resulting ex- 
pression for r is. 


- - 


2 Xi[(l+s^l-s)]^J^s) 

ik ^K2_(ik) + KQ(ik)j 


( 3 . 22 ) 


where Kq, K^_ are modified Bessel functions of the second kind. 
The subscript on T is used to indicate the Kutta condition. 
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On the other hand. If the surface has a finite thickness, 
the value of r Is diminished in a way that must be calcu- 
lated with the viscous theory. An extreme case is that 
when no vorticity is shed into the wake or f = 0. One 
may refer to this example as a case of no Kutta condition. 

It is physically the more correct result for high frequency 
finite thickness surface interactions where there is in- 
sufficient time for viscosity to establish a finite circula- 
tion. Recall (Figure 2), however, that even at low fre- 
quency there are significant departures of the total lift 
from the result produced by the Kutta condition. 

The expression for the total lift on the surface is 


L = r [ik e^^K^(ik)] -2ik j (1 - s^)^J^s)ds 

( 3 . 23 ) 


For k -► 0 the lift is proportional to V. For high fre- 
quency and no Kutta condition the lift is given by 


Lnk = -2ik 



(1 


s 2 J^s ) ds 


(3.24) 


On the other hand if the full Kutta condition is applied at 
high reduced frequency, the asymptotic form of the lift is 
given by the expression 



[d+s )/(l-s)] ^"^s) ds 



(1-s^ ) ds 


( 3 . 25 ) 


The spectral form of the total lift is of interest 
because it is proportional to the far field acoustic spectrum 
in the compact limit. Consider the compact case where the 
reduced frequency is large and the Helmholtz number < in 
(3.21)ls small or 


tc << 1 << k (3.26) 


kkl 
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The Mach number must, of course, be small for the validity 
of these relations. Then only the leading term in the series 
expression in (3.21) is Important and S(a) is given by the 
asymptotic expression 

S(a) = r ik (Tr/2k)^exp(-i7r/i|) -2ik / ^ (l - s2)^^s)ds 

y_i 

(3.27) 

With the full Kutta condition and no Kutta condition, the 
asymptotic spectrum is proportional to and respectively . 

For convective disturbances (e.g., a frozen turbulence 
pattern that interacts with the surface) the spectrum of the 
surface velocity is of the form 

^(s) = ^ (3.28) 

where is some function of the frequency. In this case 

the integral in (3.27) can be evaluated exactly and the 
acoustic spectra and lift are given by the following ex- 
pressions : 


Sk = Lk = -2 (3 25 ) 

®NK “ ^NK “ • i sln(k - tr/H) (3.30) 


and for the spectral density functions, 

|SkI^ - IL^I^ - SttI ■^1^ • 1/k (3.31) 


° ^^1^ • ~ Sln“(k - v/H) (3-32) 

k 


The Interesting feature of these results is that with a 
Kutta condition the spectrum decays monotonlcally with k 
while the no Kutta spectrum has an oscillatory decay with 
local maxima and zeroes for particular values of k as 
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shown In Figure 3* -These results suggest that a simple 
experiment could be derived to measure the acoustic spec- 
trum of a controlled convective or acoustic disturbance and 
determine by what degree the Kutta or no Kutta behavior is 
prevalent at high frequency. 

B. Vortex Fly By a Sharp Edge 

A specific example of the foregoing result is that of a 
single weak vortex convecting past a fixed surface. In that 
case the surface upwash is 


(x,t ) 


Tq X - Voot 

2tTV^ ( X - Voot ) 2 + h^ 


(3.33) 


where Fq Is the strength of the vortex and h 
vertical distance of the vortex from the plate, 
of is of the form (3*28) with 


is the 

The spectrum 



-kx 


2v^2 


where 


T = 2h/c 


(3.3^) 


(3.35) 


and k is assumed to be positive. For convecting boundary 
layer turbulence, the parameter h is of the order of the 
boundary layer thickness, the spectrum is peaked and has a 
sharp cutoff near the reduced frequency, k = 1/t, that is 
much greater than unity. Substitute '(3.34) into (3-31) and 
(3.32) and estimate the acoustic Intensity with the relation 


(2 v„/c)^ 



The final result is 


(3.36) 


2TTa3 


£2 


sin'^e 



ke”^^*^ F*dk 


(3.37) 


where 


F» = 1 


Kutta Condition 


1 


i/ I 
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= 4sln^(k - tt/H ) 


No Kutta Condition 


(3.38) 


The integration over k is easily performed with and without 
the Kutta condition. The ratio of the acoustic intensity with 
and without the Kutta condition is 


,J^K , 1 

(2 - ilT’)/(l + 


= 1/2 for T << 1 


(3.39) 


Thus the Intensity is approximately 3dB less when the Kutta 
condition is applied. With a full Kutta condition convective 
vorticity is produced in the wake with an energy level that 
is the same order of magnitude as the radiated sound. In 
the supplement to this report a full viscous problem with no 
flow is considered and a similar result is obtained. The 
vortical energy in that case, however, is confined to the 
vicinity of the edges and is not convected into the far wake 
as above. An even more drastic result was obtained in the 
recent trailing edge noise study of Howe (Ref. 13).. For a 
semi -inf inite plate the edge noise due to convected turbulence 
is typically lOdB greater when no Kutta condition is applied. 
It should be pointed out however that the semi-infinite plate 
presents a much more singular problem because of the possi- 
bility of very weak upstream decay of the surface loading. 

IV. SUPPLEMENTARY MATERIAL 

■A. Discussion of Trailing Edge Noise Models 

During the earlier part of this program an inviscld dy- 
namic model of discrete vortex fly by of a blunt (parabolic 
and rectangular) trailing edge was developed and calculations 
of the surface pressure and spectra near the trailing edge 
were made. The objective was to compare with the experimen- 
tal results of Brooks and Hodgson (Ref. 6). The Inviscid 
model accounts for two effects that are not included in the 
preceding analysis. First, the geometric shape of the 
edge (and the associated curvature of the mean flow) and 
second, the non-linear reaction of the convecting vortex to 
the surface boundary conditions. On the other hand, the 
important effect of vortex production and shedding due to 
viscosity is omitted. 
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The discrete vortex model was used in earlier work CRefs. 
1,2) to illustrate the basic coupling between the vorticity 
and acoustic modes. The model has been used by Hardin (Ref. 
14) to model the hydrodynamic origin of cavity noise. 

More recently Obermeler (Ref. 15) has- used the model to cal- 
culate the noise eiraninated due to vortex fly by of finite 
cylindrical bodies. Again, no vortex shedding is permitted 
in Obermeler 's study and noise Is calculated with the theory 
of Mohrlng (Ref. l6) that in practice differs very little 
from the theory of Powell (Ref, 17) as used by Hardin (Ref. 
li^). 


While some of the basic velocity scaling laws and fre- 
quency parameters evolve from the inviscid discrete vortex 
model, the fundamental question of the trailing edge noise 
problem is side stepped completely. That Is, how does vis- 
cosity control the vortex shedding and load alleviation near 
the edge? Howe, (Ref. 13) in a recent review of trailing 
edge noise theory has called attention to the Kutta condi- 
tion as the fundamental uncertainty in the edge noise problem 
(See Section IIIB). If one adopts a potential theory point 
of view with regard to the trailing edge fluid mechanics, 
and, furthermore, if the edge is mathematically sharp, then 
the difficulty can Indeed be reduced to an argument about the 
potential flow edge singularity or the Kutta condition. 
However, it is this authors point of view that "potential 
theory" as a fluid mechanical model is at its very worst near 
the trailing edge. Furthermore, with the slightest amount of 
edge blunting the edge singularity is eliminated and the 
Kutta condition becomes meaningless as a local vortex shedding 
criterion. More important is the fact that singularity 
arguments avoid rather than illuminate the basic physics of 
the problem. The main thesis of the present work is that 
viscosity is the missing piece of the physics and must be 
accounted for in the fluid mechanical calculation near the 
trailing edge. 

It was the realization of the importance of direct vis- 
cous effects in this study and in Ref. 8 that led to a redi- 
rection of the present effort from potential flow non-linear 
modeling to viscous flow linear modeling. Ideally, it would 
be desirable to calculate the complete nonlinear viscous fluid 
mechanics but that, of course, is an extremely difficult prob- 
lem at the present time. The present state of development of 
the linear viscous model is the following: The theory for 

the zero thickness two-dimensional surface is complete and 
some calculations have been made for low reduced frequency. 

The results are more directly applicable to unsteady airfoil 
work and are presented in Ref. 8. The evaluation of the vis- 
cous kernel function and the load calculations for high fre- 
quency is not trivial but considerable progress has been made. 
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Based on the flat plate results for low frequency it was 
found, however, that direct viscous effects at, the trailing 
edge are so strong that the calculations converge to the 
inviscid results with Kutta condition for Reynolds number 
(based on the chord) greater than 500. 

The foregoing result means that when direct viscous 
effects are included in the linearized problem, there must 
not be any geometric length introduced that forces an in- 
finitesimal or zero Reynolds number as for example at a 
mathematically sharp trailing edge. The correct formulation 
of the linearized viscous problem must include the finite 
trailing edge geometry, as discussed in Ref. 8. With a 
finite edge geometry the evaluation of the high frequency 
viscous function is even more difficult. To date only low 
frequency results have been obtain-ed for oscillating rigid 
airfoils. The load calculations agree more closely with ex- 
perimental results (Ref. 10) (both in magnitude and phase) 
than the flat plate calculations with Kutta condition. Based 
on the low frequency results it is concluded that even 
greater departures from the flat plate theory will result as 
the frequency is increased. The fundamental linearized 
problem, then, is to solve for the viscous load distribution 
on a surface with finite geometry for frequencies well into 
the acoustic regime. The basic theory is available and it 
remains a computational problem to carry out this program. 

B. Summary of the Supplemental Report 

In Ref. l8, explicit viscous calculations have been made 
for an oscillating airfoil problem. The results are extremely 
important both for the conceptual understanding of noise- 
surface interaction problems and for the quantitative cal- 
culation of interactive noise. The problem was motivated 
by the experiment of Brooks (Ref. 5). An airfoil section 
was vibrated at Helmholtz numbers of 1.3 and 2.04 in an 
anechoic chamber. The surface pressure and motion and the 
radiated pure tone sound were measured simultaneously. The 
surface measurements were also used in the "inviscid" Kirchhoff 
Integral relation to calculate the radiated sound. The calcu- 
lation overestimated the sound by 2 to 5 dB depending on 
direction. Also the measured chordwlse surface pressure dis- 
tribution has a flat behavior near the edges while aerodynamic 
theory predicts a square root decay of the edge loading. A 
more detailed description of the experiment may be found in 
(Ref. 5). 

The theoretical work in the supplementary report has 
shown that the measured surface load and sound field are 
mutually consistent if the effect of viscous dissipation at 
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the edges Is taken into account. The Invlscid Kirchhoff in- 
tegral has acoustic energy conservation built in so that all 
of the work done by the vibrating airfoil is propagated to 
the far field in the form of sound. When the same surface 
load is used in a viscous incompressible calculation the sur- 
face work appears in the near field (close to the edges) as 
kinetic energy associated with the vorticlty formation due 
to viscosity. For a Helmholtz number of order unity and a 
near quadrupole distribution of the surface load it is shown 
that the viscous energy dissipation is nearly equal to that 
propagated as sound. 


For Helmholtz numbers less than unity (into the compact 
regime) it is expected that even more energy will be trans- 
ferred into the vortlcity mode. Invlscid acoustic theory 
in the form of Kirchhoff integrals will in general overpre- 
dlct the noise in the compact regime. In each application 
an estimate of the dissipation should be made with an incom- 
pressible viscous analysis like that in the supplementary 
report. A correction to the invlscid Kirchhoff calculation 
can then be made. For very small Helmholtz numbers it would 
appear that there is no alternative to a complete compressible 
viscous calculation. Furthermore, when there is a finite 
convection velocity the surface work transferred into vorti- 
city is greatly enhanced and it is even more important to 
understand and calculate the effect of viscosity. 

V. CONCLUSIONS 

The most Important conclusion of this study is that vis- 
cosity is an important missing physical element in aeroacous- 
tic theories where Interactions with a surface with edges is 
involved. Even in the absence of a mean flow there can be 
an appreciable viscous effect near the edge with an associ- 
ated loss of energy to the acoustic mode. With a finite 
convection velocity, the magnitude of the circulation around 
the surface is determined uniquely by viscosity. The use of 
the Kutta condition will in general overestimate the circula- 
tion which means that more vortlcity is shed into the wake and 
less noise is generated. There is no aprlori reason why a 
potential flow singularity cannot occur at the trailing edge. 

It is much more important to estimate the magnitude of the 
load distribution correctly over the entire surface and this 
can be accomplished with the viscous theory. It is important 
that finite geometry and boundary layer displacement thick- 
ness effects be accounted for in the viscous calculation. 

The acoustic theory developed in this report and the associated 
report on viscous thin airfoil theory (Ref. 8) is the correct 
approach for carrying out the surface load calculation. The 
difficult practical step is the unsteady aerodynamic load cal- 
culation for high reduced frequency - a problem that has per- 
plexed the unsteady aerodynamlsclst for many years. 
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Considerable progress has been made during the course of this 
program but much work remains to be done before routine cal- 
culations can be made. 
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